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We propose an alternative mechanism of baryogenesis in which a scalar baryon undergoes a 
percolating first-order phase transition in the early Universe. The potential barrier that divides the 
phases contains explicit B and CP violation and the corresponding instanton that mediates decay 
is therefore asymmetric. The nucleation and growth of these asymmetric bubbles dynamically 
generates baryons, which thermalize after percolation; bubble collision dynamics can also add to 
the asymmetry yield. We present an explicit toy model that undergoes bubble baryogenesis, and 
numerically study the evolution of the baryon asymmetry through bubble nucleation and growth, 
bubble collisions, and washout. We discuss more realistic constructions, in which the scalar baryon 
and its potential arise amongst the color-breaking minima of the MSSM, or in the supersymmetric 
neutrino seesaw mechanism. Phenomenological consequences, such as gravitational waves, and 
possible applications to asymmetric dark-matter generation are also discussed. 



I. INTRODUCTION 

The standard model is incomplete: it does not accom- 
modate the observed baryon asymmetry and therefore 
new physics is required. Substantial effort has been de- 
voted to constructing theories that dynamically generate 
this asymmetry and some prominent contenders include 
GUT baryogenesis [T] , electroweak baryogenesis [2] , ther- 
mal leptogenesis [3] , and Affleck-Dine baryogenesis [U [5] . 
This paper proposes a new mechanism, which we dub 
'bubble baryogenesis'. 

Like the Affleck-Dine mechanism, our setup employs a 
complex scalar baryon <f>, represented in a polar decom- 
position as 



<t>(x) = R(x)e l0(x \ 



(1) 



where R(x) and 6{x) are four-dimensional real scalar 
fields. Under baryon-number transformations U(1)b, 4> 
rephases, R is invariant, and shifts. The charge density 
of <p is identified with the number density of baryons 



ub = Im((/>* 



R 2 0, 



(2) 



so a baryon asymmetry is present in field configura- 
tions that have 'angular momentum' in field space. Con- 
straints on B violation today imply that <p is currently 
at the origin of field space — so as not to spontaneously 
break B — and that the potential there has approximate 
U(1)b — so as not to explicitly break it. In the early 
Universe, however, we take <f> to be displaced from this 
minimum, to a place in the potential where B violation 
is more substantial. The observed baryon asymmetry is 
dynamically generated during the field's journey towards 
the origin. 

In the Affleck-Dine mechanism <fi evolves classically, 
relaxing uniformly towards the £>-symmetric minimum. 
B-violating potential terms torque 4> during its evolution, 
and so instead of moving in a straight line through field 
space, 4> takes a curved trajectory; 4> develops non-zero 
and consequently non-zero B. The phase transition from 
the i?-violating vacuum in the past to the £?-symmetric 



vacuum today is second-order or higher-order, and the 
end result is a spatially homogeneous condensate carrying 
a non-zero baryon asymmetry. 

But what if there is no classical trajectory connect- 
ing <f> to the symmetric minimum? Bubble baryoge- 
nesis occurs when <p evolves via bubble nucleation — 
either through quantum tunneling or thermal excitation. 
Spherical bubbles of true, i?-symmetric vacuum nucle- 
ate inside the false B-violating background. The bub- 
bles expand, collide, and eventually percolate; the phase 
transition completes when the entire Universe is in the 
i3-symmetric phase. During this process, baryons are 
produced through two distinct mechanisms. First, just 
as cf> receives a torque in Affleck-Dine, the instanton that 
mediates bubble nucleation also receives a torque from 
B-violating interactions. Consequently, the bubble wall 
takes a curved trajectory through field space, and it 
therefore accumulates B as it expands. Second, when the 
bubble walls collide, <j> can be excited back into a region 
of the potential where B- violating terms are large, gener- 
ating additional baryon asymmetry. In bubble baryogen- 
esis, the phase transition is first-order, and the end result 
is a spatially inhomogeneous distribution of baryons. Af- 
ter percolation, the baryon asymmetry is assimilated into 
the thermal plasma of the early Universe. 

Like bubble baryogenesis, electroweak baryogenesis 
also relies on a first-order phase transition in the early 
Universe. In that case, however, the tunneling scalar is 
the Higgs field, and baryon number is generated indi- 
rectly through scattering off bubble walls. 

In Sec. [TlJ we outline the basic elements of bubble 
baryogenesis and present a general analysis of the vac- 
uum structure, nucleation rate, asymmetry generation, 
bubble collisions, and washout. In Sec. IIIIl we define 



an explicit toy model, outline its cosmological history, 
and ascertain the final baryon asymmetry. Some more 
realistic examples — involving the neutrino seesaw mech- 
anism, and color breaking minima — are then presented 
in Sec. |IV| We discuss phenomenological signatures in 
Sec.lV] and conclude in Sec. IVII 
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FIG. 1: A potential that evolves with time so as to give a perco- 
lating first-order phase transition; the radial part of our potential 
has this behavior. We define a 'clock' parameter z that indexes 
the shape of the potential. At early times, z > 1 and the true 
minimum is B-violating and at <j> ^ 0. At late times, z < and 
the only minimum is B-symmetric and at rf> = 0. Vacuum decay is 
possible in the interval 1 > z > 0. 



II. GENERAL CONSIDERATIONS 

Of the potential for the scalar baryon V (</>), we require 
two features: first, that it vary with time in such a way as 
to yield a percolating first-order phase transition in the 
early Universe; and second, that the potential accommo- 
date explicit B- and CP-breaking dynamics. We discuss 
how these two criteria can be accommodated in Sec. Ill Al 
and Sec. |IIB[ respectively. Afterwards, we tackle the dy- 
namics of baryon production, which occurs at two times: 
first, at nucleation, which we discuss in Sec. II C[ and sec- 
ond, at collision, which we discuss in Sec. |IID[ ~Once gen- 
erated, the baryon asymmetry must persist, and migrate 
into the standard-model sector; we discuss the washout 
and decay of particles in Sec. |II E| 



A. Vacuum Structure and Tunneling 

A potential that achieves a first-order phase transition 
must evolve with time as in Fig. [I] In the early Universe, 
the stable minimum is i?-violating and at <f> ^ 0. As time 
evolves, the energy density of this i?-breaking vacuum 
increases until it is no longer the true vacuum; the lowest 
energy vacuum is now _B-symmetric and at (j) = 0, and 
bubble nucleation begins. At even later times, the energy 
density in the B-breaking vacuum increases so much that 
the minimum disappears entirely, ensuring that no region 
of the Universe remains stuck there. For the transition 
to be first-order, the bubbles must percolate before this 
time. 



It is convenient to introduce a dimensionless 'clock' pa- 
rameter z that characterizes the evolution of the poten- 
tial. The two most important events — when the minima 
become degenerate and when the i?-breaking minimum 
disappears — are taken to be at z — 1 and z = 0, respec- 
tively, as in Fig. [I] So, 



z > 1 stable vacuum at (f> ^ 
1 > z > metastable vacuum at <f> ^ 
> z no vacuum at <fi =/= 0. 



(3) 



In bubble baryogenesis, z decreases monotonically; it is 
a reparameterization of time. 

Bubble nucleation occurs in the window 1 > z > 0. 
Tunneling is mediated by an instanton, a solution to 
the Euclidean equations of motion with a single nega- 
tive mode. The rate per unit volume per unit time is 
given by 



r(z) = K(z)e 



-AS(z) 



(4) 



where A" is a determinant factor of order the character- 
istic mass scale of the potential to the fourth power, and 
AS* is the difference in Euclidean action between the in- 
stanton and the false vacuum [B], At z = 1, tunneling is 
forbidden, AS is divergent, and T is zero. At z = 0, AS 
is zero and semiclassical tunneling through the barrier is 
overcome by classical evolution down the potential. 

After nucleation, the energy difference across the bub- 
ble wall causes it to accelerate outward, rapidly ap- 
proaching the speed of light. Bubbles nucleate, expand, 
and collide until — at percolation — the entire Universe is 
in the true vacuum. When does percolation occur? 

Pick a point in space at time z. The expected number 
of bubbles N that have overlapped this point is 



N(z) 



V(z,z')T(z')^dz' 



(5) 



where V(z, z') is the three-volume at time z' of the past 
lightcone that emerged from our point in space at time z, 
and dt/dz' is a Jacobian factor. The integrand therefore 
is the probability that a bubble nucleates at a time z 
in the right position to convert our point to the true 
vacuum, and it is integrated over all past z. Percolation 
occurs at a time z„ satisfying 



1~ N(z*), 



(6) 



so that at least one bubble has nucleated in the past 
lightcone of every point in space. 

As long as z is changing slowly enough with time, then 
N(z*) ~ T*H~ 4 , where throughout the * subscript will 
denote a quantity evaluated at z — z*. The integral is 
dominated by values of z close to z* and the percolation 
condition can be approximated by 



(7) 
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FIG. 2: Depiction of a simplified bubble collision geometry. At percolation, there is typically one bubble per Hubble volume and the 
bubbles nucleated one Hubble time earlier. Since the surface area of a bubble scales as H* 2 and the volume scales as H* 3 , the baryon 
number density scales as H* . 



This means that bubbles typically nucleate one Hub- 
ble time before percolation, with roughly one bubble 
per Hubble volume at percolation, as shown in Fig. [2] 
Though some bubbles do nucleate before this time, the 
rate is too small to induce percolation. 

A time- varying potential of the form in Fig. [I] can arise 
naturally in two ways, depending on whether percolation 
completes before or after reheating. If the phase tran- 
sition occurs before reheating, then a direct coupling of 
the scalar baryon cf> to the inflaton field will give rise to a 
time-dependent effective potential. This is the same type 
of coupling that is used to generate the evolving potential 
in the Affleck-Dine mechanism [7j. From an effective- 
field-theory standpoint, such couplings are mandatory 
unless forbidden by a symmetry, and while they are of- 
ten non-renormalizable they can nonetheless play an es- 
sential role in the physics. In this scenario, the phase 
transition takes place between the end of inflation — so 
as not to dilute the baryons — and reheating. During this 
matter-dominated period the inflaton is oscillating about 
its minimum, but has not yet decayed to standard model 
particles. 

If the phase transition occurs after reheating, then a 
direct coupling of cf> to the big bang plasma will give rise 
to a time-dependent thermal correction to the effective 
potential. The same couplings that allow (j> to decay — so 
that the baryon asymmetry can migrate to the standard 
model sector — can generate such terms. Note that simi- 
lar thermal effects give rise to the first-order phase tran- 
sition in electroweak baryogenesis [5], with the crucial 
difference that the relevant scalar field there, the Higgs 
boson, is not charged under U(l) b- In electroweak baryo- 
genesis the purpose of the first-order phase transition is 
merely to provide an out of equilibrium environment for 
particle and anti-particle scattering processes. 

While bubble baryogenesis can occur in either scenario, 
the models we study in Sec. |III| and Sec. |IV| are in the 
former category, where it is easier to suppress thermal 
washout. 



B. B and CP Violation 

The Sakharov conditions [3] state that successful 
baryogenesis requires both B- and CP- violating dynam- 
ics. Under B and CP transformations the angular field 
component transforms as 

-A + const (8) 
^ -0. (9) 

By Eq. |8]) , B violation requires a potential that violates 
the shift symmetry on 0, i.e. carries explicit dependence 
on 0. Such terms are necessary for asymmetry genera- 
tion because in their absence the field has no reason to 
move in the ^-direction of field space, so by Eq. |2| no 
asymmetry is produced. In the Affleck-Dine mechanism, 
these P-violating terms torque the field in the (9-direction 
on its journey back to the origin; in bubble baryogene- 
sis, they force the instanton, which solves the Euclidean 
equations of motion, to arc in the ^-direction as a func- 
tion of spacetime. 

By Eq. jS), CP- violation requires cither potential cou- 
plings with complex phases or spontaneous breaking by 
an initial cj> localized at a CP-odd minimum. CP viola- 
tion is necessary for asymmetry generation because in its 
absence, though the potential can exert a torque, <f>'s tra- 
jectory is just as likely to curve in the -l-tf-direction as it is 
to curve in the — ^-direction. That is, in a CP-conserving 
theory, two CP-conjugate instantons contribute equally 
to the path integral. The percolating transition would 
therefore be comprised of an equal number of bubbles 
with positive and negative P, which average out to a 
P-symmetric Universe. 

Explicit CP violation breaks the degeneracy between 
these two CP-conjugate instantons. For example, one of 
them may disappear entirely if it is no longer a saddle 
point of the action. Alternatively, both CP-conjugate 
instantons can persist, but the one with a larger associ- 
ated Euclidean action will be exponentially subdominant 
to the process of vacuum decay. This will be true in the 
models we consider here, so we will only be concerned 



FIG. 3: Sample instanton profiles for the fields R and 6 as a 
function of radial coordinate p in a potential that breaks U(1)b- 
False vacuum values for the fields are indicated with dashed gray 
lines; the true vacuum is at R = 0. The center of the bubble (p = 0) 
is on the left. 



with the dynamics of the dominant instanton contribu- 
tion. 

In general, it is useful to characterize the degree of 
B- and CP-violating effects with a dimensionless 'effi- 
ciency' parameter e which is proportional to B- and CP- 
violating parameters in such a way that oc e. From 
an effective-field-theory perspective, e <C 1 is technically 
natural, but e ~ 1 is also allowed. 

Bubble baryogenesis generates baryon asymmetry in 
two ways. First, the instanton itself is asymmetric, which 
manifests itself as a surface density of baryons on the 
bubble walls. Second, bubble collisions excite the field 
back into the P-violating region of the potential. The 
net number density of baryons is given by a sum 

llB = "'Banstanton + collision- (10) 

We will discuss the two contributions in detail in Sec. Ill Cl 



and Sec. II D respectively, and show that, for a broad 
class of models, both of these contributions scale as 



n B 

.instantc 



collision 



(ii) 



where e is the dimensionless measure of B and CP vio- 
lation, R-p = |</>f| in the false vacuum, and P* is Hubble 
at the time of percolation. 



C. Asymmetry Generation: Instanton 

In the presence of B- and CP- violating potential 
terms, the instanton will arc in the ^-direction. We are 



FIG. 4: The trajectory (red) in field space taken by the Euclidean 
instanton solution of Fig. [3] overlayed on equipotential contours 
(purple) of the sample f/(l)fl-violating potential. To guide the 
eye, we have drawn a straight dashed line connecting the true and 
false vacua. The trajectory taken by the field is not straight; it 
curves in the S-direction. 



interested in computing both the net torque, which fixes 
the baryon asymmetry in the walls, and the bubble nucle- 
ation rate, which sets the percolation time z*. For both 
reasons, we need to find the instanton, since it charac- 
terizes the most likely bubble configuration to nucleate, 
and gives the rate via Eq. Q. Assuming ,50(4) sym- 
metry of the instanton, then the field components R(p) 
and 9{p) are functions of the Euclidean radial variable p 
alone, and the equations of motion for the instanton are 



R" + -R' - R6' 2 
P 



3 -t 
P 



2R 2 



d e V. 



(12) 
(13) 



Boundary conditions are regularity at the origin, so 
P'(0) = 0'(O) = 0, and that far from the bubble the fields 
settle into their false vacuum values, so R(oo) = Pf and 
9(oo) = 9p, where <j)p = Rpe 101 ' . Here, we are assuming 
that bubble nucleation happens by quantum tunneling 
through the potential barrier, as would be the case if 
percolation occurs before reheating. If instead, bubble 
nucleation occurs primarily by thermal activation over 
the potential barrier, then the Euclidean time coordinate 
is periodic, the 50(4) symmetry becomes SO(3) x U(l), 
and the equations of motion change correspondingly [TU] . 

The field value at the center of the bubble is near, but 
not exactly at, the true vacuum. Solutions are found by 
adjusting the field value at p = so that the bound- 
ary conditions are satisfied at p — s- oo; that is, we apply 
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FIG. 5: The total B that results from a numeric simulation of an 
expanding 1+1-dimensional domain wall. At nucleation, the bub- 
ble wall has B = 0; but as the wall accelerates, B rapidly asymp- 
totes to [i 2 , indicated by the dashed gray line. In higher dimensions 
than 1+1, total B would grow with the surface area of the bubble, 
but in 1+1-dimensions, 'surface area' is constant. The simulation 
was run in an expanding FRW background; gravitational expansion 
does not affect the value of fi 2 . 



Coleman's overshoot/undershoot algorithm generalized 
to two scalar field directions. A sample instanton, and 
its curved trajectory through field space, are shown in 
Figs. |] and U 

To estimate the extent of the curving, consider a simple 
potential set by only two parameters: m, the character- 
istic mass scale of the potential in the indirection, and e, 
a dimensionless parameter that characterizes the degree 
of the S-violation. Parametrically, 



d e V 
V 



(14) 



Because m is the only dimensional parameter, i?p ~ m, 
and the instanton solution varies in p o n sc ales of order 
p ~ m _1 . The parametric scaling of Eq. ( 13 ) is then such 
that 



6 1 [e ~ m ) ~ em 
AO = 0(oo) - 0(0) ~ e. 



(15) 
(16) 



To determine the 0(1) factors here requires finding the 
instantons numerically as above. 

Given instanton profiles R(p) and 9(p), the evolution 
of the bubble post nucleation follows from analytic con- 
tinuation of the instanton from Euclidean to Minkowski 
signature. The typical size of a bubble at the time of 
nucleation is much smaller than Hubble, so we can ig- 
nore the expansion of the Universe and simply continue 
p — > yjr 2 — t 2 , where r is the radial coordinate away from 
the center of the bubble, and t is time. At nucleation, 
t = and so p — r, and the field profile that nucleates is 
a slice through the center of the instanton. 

Because the bubble nucleates at rest, 9 = B = 0. How- 
ever, as the wall accelerates outwards, spacetime points 
in the wall traverse an angle in field space A9 in less and 
less time, so the baryon density inside the wall grows 



and grows. At the same time, as the bubble expands, 
the thickness of the wall becomes Lorentz contracted, so 
the baryon density is supported on a smaller and smaller 
region. As we will now show, these two effects cancel at 
late times, and the accelerating bubble wall asymptotes 
to a constant surface density p 2 of baryons. 

To compute the baryon asymmetry contained in a sin- 
gle bubble wall, we integrate Eq. ([2| on a fixed time slice 
t = t long after nucleation 



B, 



instanton 



R 2 9d A x. 



(17) 



where we plug in the analytically continued classical in- 
stanton profile R(p) and 9(p); we are working in the semi- 
classical approximation where loop corrections to this for- 
mula are small. Using spherical symmetry, this becomes 



B 



instanton 



= 47T 



4tt 



R 



(Vr 2 - r 2 ) 2 9 (Vr 2 - r 2 ) 



l dr 



-R(p) 2 9'(p)ry/p^dp, (18) 



where the prime indicates a derivative with respect to p. 
In the first line, the integral is taken from rtooo because 
analytic continuation of the instanton only gives the field 
profile outside of the light-cone; inside the light-cone, the 
field relaxes towards the B-symmctric minimum, produc- 
ing negligible baryons. The second line is obtained by 
changing integration variable from r to p = \/r 2 — t 2 . 
Lastly, we add the approximation that r is a long time 
after nucleation, considerably bigger than the size of the 
bubble at nucleation. Since R 2 9' dies off exponentially 
at large p, r 3> p over the region where the integrand has 
support, and. 



D 



instanton 



47TT 2 X 



-R(pf9'(p)dp, (19) 



which is the bubble surface area at late times (47rr 2 ) 
multiplied by the number of baryons per surface area 



-R(p) 2 9'(p)dp 



R 2 d9, 



(20) 



Because it is a line integral along the instanton field tra- 
jectory, p 2 is constant in time. In the spirit of Eq. (151, 
a parametric estimate is that 

P 2 ~ eRl (21) 

Note that the sign of pi 2 depends on the direction in which 
9 arcs, which in turn depends on the imaginary phases 
in the potential. 

Fig. [5] shows the total baryon number of a 1+1- 
dimensional expanding bubble. The stationary bubble 
wall carries zero baryon asymmetry, but as it accelerates 
up to the speed of light the number of baryons per surface 
area of the bubble wall grows and quickly asymptotes to 
p 2 . The integrated number of baryons scales with the 
surface area; in 1+1 dimensions the 'wall' is point-like, 
and its 'surface area' is constant. 
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FIG. 6: Sample field profiles of a collision event at three different 
times. The gray dashed lines indicate the value of ±0p. Top panel: 
two boosted domain walls approach one another. Middle panel: as 
the walls pass, they solve the free wave equation, and so linearly 
superimpose; the walls pass through each other unimpeded and the 
field between is deposited at — </>p. Bottom panel: on time scales 
of order m —1 , the field begins to respond to the potential. The 
field between the walls evolves back towards the true minimum, 
corresponding to behavior I in the text. 



Though we derived Eq. (20) in flat space, it remains 



true in an FRW background. Expansion acts globally on 
the bubble, affecting the growth of proper surface area 
with time, but it does not act locally — because the bubble 
wall is much thinner than H^ 1 , the field profile is not 
significantly affected, and so /j, 2 remains as in Eq. (20 1. 



Though the scale factor a(t) does modify the equation 
of motion, its effect is far smaller than the gradient and 
potential terms, which set the shape of the wall [TTj . 

Because the baryon asymmetry scales with the sur- 
face area of the bubble wall rather than the volume, 
the baryon number density from a single bubble dilutes 
with time, vanishing as r — > 00. The baryon asymme- 
try that is produced is carried away to infinity by the 
accelerating bubble walls. Thus, to explain the observed 
asymmetry, there must be many bubbles, and the bub- 
bles must percolate — not only to ensure that <j> is in the 
-B-symmetric minimum today, but also to preserve the 
asymmetry we have generated from escaping to infinity. 

At percolation, there is on average one bubble wall 
stretched across each Hubble volume. The expected num- 



ber of baryons per Hubble volume at percolation, there- 
fore, is the surface area of that wall (~ H~ 2 ) times 
the baryon surface density (/i 2 ). The number density 
of baryons right before collisions, therefore, scales like 



^£»,instanton 00 M 



X H* 



(22) 



In the next subsection, we will argue that these baryons 
are approximately conserved by the collision, so that 
this n^histanton contributes directly to the total tlb in 



Eq. (11) 



D. Asymmetry Generation: Collisions 

To gain insight into the dynamics of bubble collisions, 
we have run numeric simulations for a variety of different 
models. Specifically, we have studied the collisions of 
1+1-dimensional bubble walls in an expanding, matter- 
dominated, FRW background with scale factor a(t) for 
an array of different potentials. The equation of motion 
for 0, ignoring gravitational backreaction, is 



3H<f> 



a(t) 



2 = - V 'W- 



As initial conditions, we used the exact 1- 
instanton profiles. 



(23) 



1-dimcnsional 



The dynamics of Eq. ( 23 1 are complicated, but the 
moment of collision is simple. The colliding walls are 
relativistic — they are moving very fast and are very thin, 
by Lorentz contraction. This means that the time scale 
on which they cross is far smaller than both (the 
time scale on which FRW expansion acts) and m _1 (the 
time scale on which the potential acts) . Gravity and the 
potential, therefore, can both be ignored, and the field 
approximately obeys the free wave equation <f> — <f)" = 0, 
where we have rescaled x by a{t) the time of collision. 
Linear superposition of waves is an exact solution to the 
free wave equation, so the impinging walls merely pass 
through one another, and the field in between is deposited 
at — <f)p, as shown in Fig. [6j This behavior is generic: as 
long as the walls are moving fast enough, the field value 
at the intersection of the walls is — </>f, independent of 
the precise shape of the bubble wall or the structure of 
the potential. 

On longer time scales, of order m _1 , linear superposi- 
tion is no longer a good approximation and the field be- 
gins to evolve under the force of the potential. The field 
in between the crossed walls rolls down the potential and 
begins to oscillate around a local minimum. There are 
two behaviors, depending on which local minimum. 

/. Oscillation about the true minimum. The field be- 
tween the walls, deposited on the other side of the po- 
tential at —<t>Fi is no longer in a vacuum state; under the 
force of the potential, it evolves back towards the origin, 
as shown in the bottom panel of Fig. [6] By the time the 
field reaches the true minimum, it has lost enough kinetic 
energy to gradients and Hubble friction that it cannot es- 
cape; it oscillates, Hubble friction damps the oscillations, 
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FIG. 7: Re(</>) as function of space and conformal time rj during a 
bubble collision showing behavior I. Using conformal time ensures 
lightcones travel at 45 degrees. Domain walls separating the false 
and true vacua emerge from the top and bottom of the plot and 
collide at x = 0. The walls cross and the field is locally deposited 
at — <t>-p, from which it evolves back toward the true vacuum. The 
value of ~Re(rf>(x = 0)) is shown in red. Several localized field exci- 
tations related to oscillons are visible at late times. 



FIG. 8: The baryon number density as a function of space and 
conformal time 77 for the same collision as Fig. [7] Initially, the 
asymmetry is carried entirely on the walls; as they accelerate the 
baryon density at the center of the wall grows and the thickness 
of the wall shrinks. At the collision, the walls move through each 
other, roughly conserving baryon number, and the evolution of the 
field from — <pp back to the origin generates new baryons around 
x = 0. The oscillons of Fig. [7] carry non-zero B. 



and eventually the field settles into the true minimum. 
Fig.[7]shows a collision that illustrates this behavior, and 
Fig. [8] depicts the corresponding baryon number density. 

As the field at the collision site evolves back towards 
the origin, it moves through a i?-violating region of the 
potential, so a second wave of baryon generation is taking 
place. Fig. [9] shows the integrated baryon number as a 
function of time, which in d + 1-dimcnsions is 

B(t) = a{t) d J n B (t, x) d d x. (24) 

Before the collision B(t) was constant, as in Fig. [5] but it 
makes an abrupt jump upwards at the moment of colli- 
sion. The evolution of the field from — <pp generates new 
baryons inside the collision lightcone, visible in Fig. [8] 
A simple estimate can be made for the baryon number 
generated during this evolution. 

The process in which the field in the collision region 
evolves from — <p-p to the origin can be thought of as a 
localized Affleck-Dine condensate forming and dissolving 
at the collision site. If the field takes a time At to evolve 
from — 4>p to 0, then the spatial width of this condensate 
is of order At, since the bubble walls propagate at nearly 
the speed of light. The Affleck-Dine mechanism, were 
it to occur in this potential, would generate a number 
density of baryons which scales like R 2 9 ~ ei?|m. Mul- 
tiplying this by the width of the condensate and by the 
surface area of the collision site per unit volume at perco- 
lation gives the expected number density of baryons 



generated by the collision as 

»B .collision ~ eR%m x AtH* ~ eRpH*, (25) 

where At ~ 1/m since m is the characteristic scale of the 
potential. The contribution to tib from bubble collisions 
has the same parametric dependence as the contribution 
from the instanton. 

After the collision, the energy in the field and the walls 
dissipates, and the asymmetry spreads; the way in which 
this happens is potential-dependent. Potentials which 
grow more slowly than quadratically at the origin, ad- 
mit two related non-topological solitons that can tem- 
porarily trap energy and baryon number [12[ 113] . Field 
excitations that move solely in the (^-direction, locally or- 
biting the origin, are called Q-balls O [TS] ; their charge 
contributes a centrifugal term to their effective potential 
that makes them absolutely stable. Field excitation that 
move solely in the i?-direction, locally oscillating along 
a line through the origin, are called oscillons [IB]; these 
excitations are long-lived, but not eternal. Our collisions 
produce a hybrid: it oscillates predominantly in the in- 
direction, and can therefore be seen in Fig. [7J but it also 
carries non-zero baryon number, and can therefore be 
seen in Fig. [§} In other words, the field is locally exe- 
cuting very elliptical orbits about the origin. A number 
of these hybrids are visible: a stationary one emerges 
from the collision site, and several boosted ones fall off 
the wall as it propagates. Because these non-topological 
field configurations probe larger field amplitudes, they 
are still sensitive to U(1)b violation, as we discuss in the 
next subsection. 




FIG. 9: The total baryon number B as a function of conformal 
time for the same collision as Fig. [7] Superimposed, in red, is a plot 
of Re(0(ic = 0)) during the collision. B is flat before the collision, 
as in Fig. [5] At the collision, the field is deposited at — <pp and 
evolves back towards the origin; during that evolution, B surges. 



77. Pockets of a false vacuum. For certain potentials 
instead of oscillating around the true vacuum, the field at 
collision site ends up in a false vacuum. This can happen 
in two ways. First, if the potential at — 0p is very sloped, 
then the field can overshoot the true minimum and land 
back in the original false vacuum, as was noticed in the 
earliest simulations of bubble collisions [17] ■ Second, if 
the potential happens to have an additional local mini- 
mum at — 4>f, then the field never has to evolve anywhere, 
since it is already in a false vacuum [18] . The presence 
of an additional minimum at — <f>-p does not necessarily 
require tuning; potentials with approximate U(\)b have 
this feature automatically. 

The false vacuum provides a locally stable state for 
the field, around which it can execute small oscillations. 
Though locally stable, the field does not remain in the 
false vacuum forever. The walls, moving away from the 
collision, now have true vacuum on the outside and false 
vacuum on the inside. This induces a pressure that 
pushes the walls back towards the collision site, so the 
walls eventually slow, turn around, and re-cross on a time 
scale of order H* 1 , which is far longer than the time scale 
of oscillations about the minimum. The formation and 
collapse of long-lived pockets of false vacuum is shown in 
Fig. [lO] While the field lingers in the false vacuum, B 
is not conserved, and the fact that the field is oscillat- 
ing around the false vacuum can yield large fluctuations 
in the baryon asymmetry. In this case, the asymmetry 
is presumably still non-zero, but it is difficult to get an 
analytic handle on it, and numerical simulations are re- 
quired. 

For potentials with broken U(1)b, there is no reason 
for — </>p to be a local minimum, and the steepness of the 
potential at — (\>p determines which behavior occurs. For 
potentials with approximate t/(l)s, the physics depends 
on the nature of the false vacuum. When — <fip has a 
strong basin of attraction, which happens at larger z*, 



Re(cp) 




FIG. 10: Re(0) as function of space and conformal time r\ during 
a bubble collision showing behavior II. As in Fig. [7] domain walls 
separating the false and true vacua emerge from the top and bottom 
of the plot and collide at x = 0. At the collision, the field is 
deposited at — 0p, which is another false vacuum of the theory. 
The field oscillates about this vacuum, while the bubble walls move 
outward, slow, turn around and re-cross. This process can occur 
several times. 



collisions tend to exhibit behavior II; when —(f>F has a 
weak basin of attraction, which happens at smaller z*, 
collisions tend to exhibit behavior I. In the toy model of 
Sec. |III[ percolation tends to happen at smaller z„. 



E. Washout and Decay 

After percolation, a baryon asymmetry density of order 
~ eRpH* is inhomogeneously distributed throughout the 
Universe in the form of the <f> field. In order to explain 
observation, this asymmetry must persist and it must 
migrate to the standard-model sector. 

For the asymmetry to persist, we must avoid two types 
of washout: classical and thermal. Classical washout 
refers to depletion of the asymmetry by evolution under 
the classical equations of motion from B-violating opera- 
tors present in the early Universe. The dynamics of clas- 
sical washout depend on the dimensionality of these oper- 
ators. In the case of higher-dimension operators, washout 
can only be effective at large field values. But, since the 
expansion of the Universe and the growth of the bubble 
both tend to damp field excitations toward the origin, 
classical washout from higher-dimension interactions is 
typically evaded, as is true in the Affleck-Dine mecha- 
nism. In the case of marginal or super-renormalizable 
operators, classical washout may be active even as the 
fields damp to origin. Consider, for instance, the case in 
which the potential for </> near the origin is of the form 
m 2 (|0| 2 + e<fi 2 + h.c), where e is a dimensionless mea- 
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sure of P violation. The e term induces an 'ellipticity' to 
the potential that splits the mass eigenstates, causing the 
field to precess as it orbits the origin. The total baryon 
asymmetry, therefore, oscillates around its initial value. 
The frequency of these oscillations scales like ~ (em) -1 . 
Whether the baryon number is spread out evenly, or lo- 
calized in non-topological solitons, the effect of classical 
washout is to make it oscillate. 

For small e, this precession frequency is far lower than 
the characteristic oscillation time of the field m _1 ; the 
field must sit at the origin through many oscillations if 
the asymmetry is to appreciably change. As long as the 
<f> condensate decays before this time, the asymmetry is 
preserved. Besides, even if it does not decay in time, 
the asymmetry oscillates around its initial value, it does 
not damp; unless there is a conspiracy between the decay 
time and the oscillation time, the final asymmetry will be 
an 0(1) fraction of the initial asymmetry. Alternatively, 
in certain models P-violation is 0(1), in which case nu- 
merical simulation is necessary to evaluate the degree of 
washout. 

Thermal washout refers to depletion of P that occurs 
after the asymmetry is absorbed into the plasma of the 
early Universe, through scattering processes which in- 
volve P-violating interactions. Such scattering tends to 
restore chemical equilibrium and therefore deplete the 
asymmetry. Even if these interactions arise from higher- 
dimension operators, they can still be significant since 
the associated scattering rates grow with temperature. 
However, as we will discuss in Sec. our models evade 
thermal washout because P-violation is sourced by inter- 
actions between (f> and the inflaton. P-violation occurs 
when the temperature is small, before the inflaton has 
decayed; when the inflaton decays and the temperature 
becomes large, the P-violating interactions have shut off 
and the asymmetry is frozen in. Thermal washout is 
avoided because P-violating interactions and the ther- 
mal plasma are never present at the same time. 

Finally, for the baryon asymmetry to migrate to the 
standard-model sector, the <f> condensate must decay. 
A direct coupling of <f> to some operator comprised of 
standard-model fields suffices to transfer the asymmetry; 
the decay rate of 4> depends on the strength of this cou- 
pling. For homogenous fields, the decay rate of a conden- 
sate is more or less the same as the decay rate of <j) quanta 
in the vacuum [19] . In bubble baryogenesis, however, the 
field configuration is highly inhomogeneous. In partic- 
ular, because the bubble walls are boosted to near the 
speed of light shortly after nucleation, one might worry 
that these field fluctuations will be long-lived due to a 
Lorentz boost factor. However, as discussed earlier, af- 
ter the collision these boosted field configurations are not 
solutions to the equations of motion so they broaden and 
dissolve. As the quanta become softer, they can decay. 
In certain models of bubble baryogenesis, the decay of <j> 
can be very fast — much faster than the Hubble parameter 
at the time of percolation. In such cases, the <fi conden- 
sate decays to particles nearly instantaneously after the 



nucleation event. Afterwards, there is no classical field, 
meaning classical washout is straightforwardly evaded. 

III. TOY MODEL 

Our discussion in Sec. [TT] was framed rather broadly, so 
in this section we now study a concrete setup. A work- 
ing model of bubble baryogenesis must accommodate the 
following criteria: 

• Percolation. The Universe must efficiently transit 
from the false to true vacuum. 

• Asymmetry. The theory parameters must accom- 
modate the observed baryon asymmetry. 

• Pcrturbativity. The effective couplings cannot blow 
up and all energy scales are bounded by the cutoff. 

• No Washout. B- violating effects, classical or ther- 
mal, must be under control. 

In this section we present a toy model which satisfies 
these criteria. As we will see, despite its simplicity the 
toy model may actually be phenomcnologically viable. 

A. Model Definition 

Our toy model is defined by a potential of the form 

V(R,9) = VoW+eVtiW), (26) 

where e is a small parameter characterizing B- and CP- 
violating interactions. The P-symmetric part of the po- 
tential is 

V = m 2 H 2 -A|0| 3 + AH 4 . (27) 

In general, the shape of the potential will vary in time 
due to couplings between <fi and the inflaton x> whose 
vacuum expectation value is time-dependent. This oc- 
curs in Affleck-Dine models of baryogenesis, where such 
couplings induce Hubble-dependent parameters in the ac- 
tion. For simplicity, we take A and A to be constant, but 

m 2 = rh 2 -^, (28) 

where p = 3P 2 mpj is the energy density of the Universe. 
Here we require A 2 > and m 2 > so that P is sponta- 
neously broken at early times but restored in the present 
day. In a supersymmetric context, such a p dependence 
would originate from |x| 2 |(/>| 2 /A 2 in the Kahler potential. 

If the couplings between </> and x are B~ an d CP- 
violating, then V\ will also contain time-dependent 
terms. At least two P-breaking operators are required — 
otherwise all CP phases can be removed by a field re- 
definition. As discussed in Sec. |II E| if P violation is 
higher-dimensional, e.g. (f> 5 and </> H , then explicit break- 
ing is localized far from the origin and classical washout is 
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FIG. 11: AS(z) for the potential in Eq. \27\ , computed numer- 
ically (blue solid) and analytically in the thin-wall approximation 
(red dashed). 



ameliorated by Hubble damping of the field. In contrast, 
renormalizable P-violating operators, e.g. 4> 2 and (j> 3 , typ- 
ically mediate classical washout with model-dependent 
effects. Importantly, we assume that P-violating inter- 
actions are sourced by the inflaton alone, so P is re- 
stored after reheating — thus, thermal washout is evaded. 
This setup can be easily engineered by invoking addi- 
tional symmetries under which both (f> and \ transform. 

Let us outline the cosmological history of this model. 
We begin during inflation, when H and therefore all 
model parameters are constant in time. During this 
epoch, to 2 < and the field resides at <fi ^ 0. After 
inflation ends, the inflaton starts to oscillate around its 
minimum and the Universe shifts from vacuum-energy 
domination to matter domination. As H decreases, m 2 
eventually becomes positive, growing monotonically un- 
til the 0^0 vacuum becomes metastable. Once the 
nucleation rate rises sufficiently, percolation occurs. The 
P- and CP-violating interactions in the potential cause 
the nucleated bubble walls to accumulate baryons; sub- 
sequent bubble collisions yield an additional asymmetry 
component. 

In the subsequent sections, we will analyze the vacuum 
structure, nucleation rate, and asymmetry generation in 
this model. 



B. The Instanton 

Neglecting effects proportional to e, the Euclidean ac- 
tion for this theory is 



S = 2tH 



p 3 dp 



dR 



m 2 R 2 



AR 3 + X 2 R 4 



(29) 



While the model parameters vary in time, they do so on 
scales H <C 1/m, so we treat them as constant in our 
analysis of the instanton. It is convenient to transform 



s = 


S/X 2 


R = 


Rm/X 


A = 


AraX 


P = 


p/m, 



where 



S = 2ir 2 / p 3 dp 




R 2 - AR 3 + R 4 



(30) 



(31) 



As a result of our change of variables, the rescaled Eu- 
clidean action S is a function of A alone. 

As discussed in Sec. |II A[ the variable z is a convenient 
reparametrization of time in which 1 > z > corre- 
sponds to the epoch in which tunneling is allowed. We 
can express A in terms of z as 



A(z) = V32/9 + z(2- a/3279), 



(32) 



where 1 > z > maps onto the range 2 > A > a/32/9. 
Within this interval, there is a true and false vacuum 
located at 



i?T — , i?p 



3A 

IT 



1 



1 



32 \ 
9A 2 ) 



(33) 



Next, we evaluate AS* by solving the associated Eu- 
clidean equation of motion for an SO (A) symmetric 
ansatz R subject to the initial condition OR/ dp = at 
p = 0. Solving for AS* numerically, we find that for 
z < 0.25, AS" is very well fitted by the function 



AS = 431.52 



0.679 



8139.4z 



2.27 



(34) 



As z — ¥ the phase transition shifts from first-order to 
second-order. As z —> 1 the bubble becomes thin-walled. 
From Fig. |11[ it is clear that our numerics agree with 
analytic expressions in this regime. 

The determinant prefactor can also be straightfor- 
wardly estimated. From 20] we can compute the K fac- 
tor in Eq. Q 

-1/2 



K 



AS 

4ir 2 



det'(-D + V"{<p{p))) 



dct(-D + V™(0F)) 



AS 2 m 4 
4rr 2 A 4 : 



(35) 



(36) 



where <f>(p) is the instanton solution and det' indicates the 
determinant with zero modes removed. In the second line 
we used that p varies on scales of order to -1 to estimate 
the determinant factors. 

Up till now we have neglected B- and CP-violating ef- 
fects proportional to e. Employing our numerical code, 
we have generated instanton profiles for the potential at 
finite e. Given these n um erical solutions, we can com- 
pute p 2 , defined in Eq. ( 20 ) as the measure of the baryon 
asymmetry in the bubble wall. As shown in Fig. |12[ at 
small e, these numerical resu lts match the simple esti- 
mate for p, 2 described in Eq. (21 ). 
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FIG. 12: Plot of the exact value for y? divided by our parametric 
estimate, tRp — the answers agree within O(l) factors. Here we 
used Vb corresponding to z = 0.2 and V\ = R 2 cos(2# + tt/4) + 
R 3 cos(30). 



C. Before Percolation 

Consider the cosmological history of this model leading 
up to percolation. As shown in Sec. |III B the vacuum 
structure of the theory depends solely on A, which varies 
in time with the energy density, p. Shortly after inflation, 
m 2 < and the potential has a single minimum at large 
field values. As the universe cools, eventually in 2 > 
and an additional local minimum forms at the origin — at 
this point A is divergent. As m 2 continues to increase, 
A monotonically decreases, and a first-order phase tran- 
sition becomes possible in the window 2 > A > y/32/9. 



Plugging p = 3iJ 2 m P1 into Eq. (|32|) and Eq. (|28|, we 



obtain expressions for the important physical quantities 
as a function of z, 



m(z) 
H{z) 



A 



XA(z) 



Ay^rh 2 — m(z) 2 
y/3mp\ 



(37) 
(38) 



where A{z) is defined in Eq. (32). Given the parametric 
dependences in Eq. ( 38 ) , it is clear that during the first- 



order phase transition, the dimcnsionful parameters, H, 
m, and A are all well below the cutoff A and so the 
effective theory description remains valid. 

Eventually, the universe cools sufficiently that Hubble 
decreases enough that there is a first-order phase transi- 
tion. Percolation occurs when Eq. ^ is satisfied, which 
is roughly when 1\ ~ Hf; in the context of our model 
this approximation is accurate to < 15%. Recall that 
1\ implicitly depends on through the model param- 
eters in Eq. (38). The criterion for percolation can be 
rewritten as 



ASjmj 
At: 2 



,-AS, 



(39) 



this equation, and so the prefactors are only logarithmi- 
cally important. Solving for S* yields 



AS*, = log 
~ 4 log 



9m P1 AS 2 



4tt 2 A 4 (1 -m 2 /m 2 )' 



/mpi 
V A 



(40) 
(41) 



where the ellipses denote logarithms of 0(1) numbers. If 
S, C 1, then the nucleation rate is very high at the onset 
of percolation, indicating that the phase transition is bor- 
dering on second order. While this is not necessarily bad 
in and of itself, in this regime <fi evolves like a slow roll 
Affleck-Dine condensate. Of course, we will be interested 
in the first-order regime, whereby tunneling is the dom- 
inant mode of the phase transition. To achieve this, we 
require a modestly sized instanton action, so A < TOpi, 
which is to say that the higher dimension operators cou- 
pling the inflaton to <j) cannot be Planck slop operators, 
and must be suppressed by a lower scale. 

In this toy model is typically small and we are in 
the weakly first-order regime. Allowing both A and A to 
also vary with time can yield any value of z+ between 
and 1. 



D. After Percolation 

Once percolation occurs, bubbles of true vacua nucle- 
ate and soon fill the volume of space. Using Eq. (21), 



which applies for e <C 1, we find that the surface density 
of baryons on the walls is 



IK 



A 2 ' 



(42) 



The total baryon number density Ub arising from the 
initial instanton plus the subsequent bubble collisions is 
given in Eq. ( 11 ). To compute the observed baryon asym- 



metry today, we need to consider the remainder of the 
cosmological history. Because the inflaton sources the 
B- and CP- violating interactions of cf>, we require that 
the decay of the inflaton, and thus reheating, occur af- 
ter the percolating phase transition. Consequently, the 
asymmetric yield at the time of reheating is given by 



Vb 



n B #r 
sr H 2 

4A 2 mpi KyJ in 2 / m 2 



(43) 
(44) 



The exponential factor essentially fixes the solution to 



where and Tr are the Hubble parameter and the 
temperature, respectively, at the time of reheating. The 
above estimates neglect the effects of classical washout, 
but as discussed in Sec. |II E| these effects are expected 
to change the asymmetry by an 0(1) fraction and are 
model-dependent . 

To accommodate the present day observed baryon 
asymmetry, we require that rjg ~ 6 x 10~ 10 . The en- 
ergy density of the Universe at reheating is bounded 
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decays so fast into particle quanta. Alternatively, fast 
decays can occur if <fi decays promptly to additional non- 
gauge-singlet particles which in turn decay to standard 
model fields. Finally, while our discussion thus far has 
been framed within the context of B, the asymmetry can 
of course be converted into L via the operators 4>LH U , 
<t>LLE, or <j)QLD, or into a dark matter asymmetry via 
similar interactions. 



IV. MORE REALISTIC MODELS 

In this section we discuss possible realizations of bub- 
ble baryogenesis in more realistic contexts. The examples 
here are supersymmetric, but we emphasize that this is 
not a requirement for bubble baryogenesis. 

A. Neutrino Seesaw 



FIG. 13: Parameter space for the toy model, fixing rjg to the 
observed baryon asymmetry today. The blue region depicts the 
allowed region after requiring that reheating occurs after percola- 
tion, ppt < p*. The red dotted lines are contours of Tr and the 
solid gray lines are contours of AS 1 *. We have fixed A = 1, e = 0.1, 
rh = A. 



from above by the energy density at percolation, so 



PR ^ P*- I n turn, this places an absolute upper 



limit on r]B from Eq. (44). Given the observed baryon 



asymmetry, this upper bound can be rephrased as a lower 
bound on rh in terms of the other fundamental parame- 
ters 



> A 



10 1U GeV x A z 

7K 



2/3 



(45) 



where for simplicity we have dropped the 0(1) factors 
from the dependence. This limit is a substantial con- 
straint on rh in this toy theory, and is depicted in the 
blue region of Fig. [13} 

Next, we briefly discuss how the asymmetric yield 
in Eq. (44) is actually transferred from <j> into stan- 
dard mode fields. In a supersymmetric context this is 
achieved, for example, by the operator UDD/M in the 
superpotential, where M is the mass scale of some con- 
nector field which has been integrated out. This is the 
lowest-dimension operator which can link standard model 
B to a gauge singlet field <f>. Given this operator, the field 
4> has a decay rate of 



r(</> -> qqq) 



1 



128tt 3 M 2 ' 



(46) 



which can be much greater than H*, the Hubble pa- 
rameter at the time of percolation. In this case the <$> 
field within each bubble of nucleated true vacuum de- 
cays very shortly after percolation. Classical washout 
from _B-violating interactions are thus minimized since 4> 



Consider the MSSM augmented by a supersymmetric 
neutrino seesaw. The superpotential is 



W = \LH U N + ^MN 2 



(47) 



where N denotes the sterile neutrinos and we have sup- 
pressed all flavor indices. Explicit violation of U(1)l is 
present in the form of the mass parameter, M. Inte- 
grating out the heavy N fields yields the active neutrino 
masses, whose sum is bounded by cosmological measure- 
ments to be less than 0.17 eV [5T]. Requiring that A < 1 
implies that M < 10 14 GeV. 

The supersymmetric F-term and D-term contributions 
to the potential are 



V F = \\H U N\ 2 + \XLN\ 2 + \XLH U + MN\ 



V D = 



g' 2 + g 2 < 



\H u \ 2 -\L\y 



(48) 



If supersymmetry breaking approximately preserves L, 
then the corresponding contributions to the potential are 
of the form 



V 



V Q + eV t 



where e is small and 



Vn = mi\L\ 



-ALH U N + h.c. 



'H, 



(49) 



(50) 



Because the parameters in V acquire contributions 
from both zero-temperature supersymmetry breaking 
and Hubble-induced supersymmetry breaking, they are 
in general time-dependent. Which specific contributions 
arise depends on the symmetry structure of the ultra- 
violet theory, which dictates the coupling between the 
MSSM fields, the inflaton, and the supersymmetry break- 
ing sector. For instance, if the inflaton carries i?-parity, 
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then L-violating terms like LH U or N 3 could be present 
inVi. 

The scalar potential is complicated, and contains a 
large number of fields and parameters. However, it is 
clear that all the required features of bubble baryogene- 
sis are present. First, the A parameter may be large in 
the early Universe and produce global minima far from 
the origin of field space. Second, these _B-breaking vacua 
are stabilized by the quartic A which can be 0(1) in order 
to achieve a sufficiently large nucleation rate. Third, L- 
and CP- violating interactions can torque the Euclidean 
instanton solution during a first-order phase transition. 
Hence, we do not expect bubble baryogenesis in this po- 
tential to differ in any qualitative way from the toy model 
in Sec. QUI 

A complete analysis of the multi-field potential of the 
supersymmetric neutrino seesaw would be non-trivial. 
However, we take note of certain simplifications to the 
theory which can reduce the potential to a solvable one. 
In particular, the lepton flavor indices serve largely to 
complicate the analysis, so the potential is greatly sim- 
plified by taking a single flavor of L and N to be the only 
fields active in the dynamics. Furthermore, the D-term 
contribution in Eq. (48) tends to fix \L\ = l-ffij 1 , which 
effectively eliminates another set of field directions. 



B. Color-Breaking Minima 

Bubble baryogenesis may also be possible within 
the context of tunneling from color-breaking minima 
[52] to the electroweak vacuum within the MSSM. If 
supersymmetry-breaking A-terms are too large, then 
deep minima can form at field values away from the ori- 
gin, inducing an instability for the electroweak vacuum. 
The vacuum dynamics are dominated by the field direc- 
tions of the top squark and the Higgses, where the su- 
perpotential is dominated by the usual top quark Yukawa 
coupling, 



W = \ 3 Q 3 H U U 3 



(51) 



Given approximately B-symmetric supersymmetry 
breaking, the corresponding terms in the potential are 
as in Eq. (|49| but with 



Vo = m 2 Q3 \Q 3 \ 2 + m 2 U3 \U 3 \ 2 +m 2 H jH u \ 2 

+A 3 Q 3 H U U 3 +h.c. (52) 

Considering the D-flat field direction |Q 3 | = \U 3 \ = \H U \, 
then the absence of color-breaking minima implies that 



A 2 <3X 2 {m 2 Q3 +m 2 U3 +m 2 H J, 



(53) 



if the electroweak vacuum is to be absolutely stable. 

For our purposes we take the opposite tack — we want 
the electroweak vacuum to be unstable in the early Uni- 
verse. Couplings to the inflaton can be easily arranged 
so that Eq. (53) fails in the early Universe, so that the 



fields reside in the true, color-breaking vacuum. As the 
Universe cools, Eq. (53) is eventually satisfied, and the 



1 Some cosmological implications of such D-flat directions are dis- 
cussed in 1231. 



fields can tunnel from the color-breaking minimum to the 
electroweak vacuum. 

Of course, to generate a baryon asymmetry, this phase 
transition requires B- and CP-violating interactions. 
The natural candidate for this is the Hubble-induced 
cubic term, UiDjDk, whose coupling can carry non- 
zero CP phases. Because this operator involves multi- 
ple squark flavors, we are then required to understand 
the tunneling from the color-breaking minimum in other 
squark directions besides the stop. We leave a proper 
analysis of this scenario for future work. 



V. OBSERVATIONAL CONSEQUENCES 

The byproduct of bubble baryogenesis is a frothy mix- 
ture of standard model baryons, inhomogeneously dis- 
tributed in the early Universe. This inhomogeneity pro- 
vides a potentially dangerous relic, since the observed 
baryon density is known to be homogeneous at the epoch 
of big bang nucleosynthesis (BBN). On average, inho- 
mogeneities are small at BBN: the percolating bubbles 
have a length scale set by H~ 1 , which is far smaller than 
(a*/aBBN)^BBN' the size of an inhomogeneity at perco- 
lation that would grow to be Hubble size at BBN. The 
potentially dangerous relic doesn't come from the aver- 
age bubble, it comes from the rare bubble that nucleates 
early enough and avoids enough collisions to grow big by 
percolation. Constraints from big bubbles were studied 
in [53] in the context of extended inflation 25j , which also 
features a first-order transition around the end of infla- 
tion. Though the context was different, the constraints 
are purely geometrical, so the analysis in [24] carries over. 
Big bubbles constrain the decay rate considerably before 
percolation to be small, so that big bubbles are exceed- 
ingly unlikely. Bubble baryogenesis is aided in avoiding 
this constraint by the fact that the nucleation rate can 
be completely shut off at early times. 

Constraints aside, the first-order phase transition also 
opens up new observational signatures, like gravitational 
waves. Bubble collisions are an efficient producer of grav- 
itational waves; numeric estimates in |26) show that as 
much as .1% of the energy released in the transition can 
end up in gravitational waves. Unfortunately, the energy 
for us is small — most of the energy density of the Uni- 
verse is in the inflaton — but the gravitational waves have 
a distinct signature which may make observation feasi- 
ble. Because the colliding bubbles at percolation have 
roughly the same size, the gravitational wave spectrum 
has a spike at [27] . This observational signature is 
distinct from Affleck-Dine. 
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Additionally, black holes might form at the collision 
sites. This is an intriguing possibility, because a cur- 
rent tension in the data of the mass distribution of early 
quasars would be alleviated by a source of primordial 
'seed' black holes [21] • Also, because these black holes 
form with a characteristic size, there could be a second 
bump in the gravitational wave spectrum from their co- 
incident evaporation. 

Lastly, 



II D bubble collisions 



as we discussed in Sec. 
can spawn non-topological solitons, like oscillons and Q- 
balls. Oscillons, though long-lived, are not typically sta- 
ble on present Hubble time-scales [2T)H3~1] ; they radiate 
energy and dissociate. Q-balls, on the other hand, can be 
stable and persist since they are charge-stabilized [321I33]. 
In Affleck-Dine baryogenesis, Q-balls typically form in 
gauge- mediated theories [34] , and may or may not form 
in gravity-mediated theories, depending on the super- 
symmetric spectrum [35, 36J. Analogous statements are 
likely true for bubble baryogenesis. If Q-balls form, their 
subsequent evolution is model-dependent: Q-balls that 
carry B will be absolutely stable if their mass-to-baryon- 
charge ratio is smaller than that of a proton; Q-balls 
that carry L, however, typically decay, since the mass- 
to-lepton-charge ratio of the neutrino is so small. 



VI. FUTURE DIRECTIONS 

Bubble baryogenesis is a novel scheme for the gener- 
ation of the cosmological baryon asymmetry. A scalar 
baryon undergoes a first-order phase transition in the 
early Universe, and a baryon asymmetry is generated by 
the process of bubble nucleation and the subsequent bub- 
ble collisions. We have presented an explicit toy model 
to illustrate the basic features of the mechanism, and in- 
troduced a handful of realistic models. In addition to 
fleshing out these more realistic theories, there exists a 
variety of interesting directions for future work. 

First, in Sec. |II C| we argued that a single bubble 
could never explain the observed baryon asymmetry be- 
cause the baryons are carried away to infinity. This is 
not strictly true: a loophole is provided in theories with 
large extra dimensions — the same loophole that boom- 
and-bust inflation [37] exploits to provide a graceful exit 
to old inflation. If a bubble nucleates smaller than the 
size of the extra dimensions, then as it grows, it wraps 
the extra dimension and collides with itself on the other 
side. The bubble wall no longer runs off to infinity; the 
self-collision preserves the baryons in the wall, and dis- 
tributes them uniformly throughout the interior of the 



bubble. 

Second, as we discussed in Sec. II A there is the pos- 
sibility of bubble baryogenesis after reheating. Ther- 
mal effects from the big-bang plasma can induce time- 
dependent couplings which give rise to a first-order phase 
transition; the thermal plasma also assists the transition 
by enhancing the nucleation rate. To achieve such a sce- 
nario requires engineering the appropriate thermal po- 
tential and avoiding thermal washout. 

Third, models considered in this paper had only a sin- 
gle instanton mediating decay, which imposed a relation- 
ship between the efficiency of the phase transition and 
the resulting asymmetry This need not be the case. 
Consider a true and false vacuum connected by two in- 
stantons: a dominant one that is purely radial so that 
it generates zero baryons, and a subdominant one that 
arcs so that it alone is responsible for generating baryons. 
Physically, this would correspond to a percolating phase 
transition in which the vast majority of nucleated bubbles 
are .B-symmetric, but some small fraction are asymmet- 
ric. The smallness of the asymmetry in such an example 
would arise not from small e, but from the exponential 
suppression of the subdominant instanton. Such a model 
may suffer from fine-tuning issues because the baryon 
asymmetry would be so sensitive to the Euclidean action 
of the subdominant instanton. 

Finally, in recent years there has been a resurgence of 
interest in so-called asymmetric dark matter, where the 
dynamics of baryogenesis and dark matter genesis are 
linked |38j . Such a linkage can arise naturally in hidden- 
sector theories in which dark matter has a U(1)dm, an d 
can have phcnomcnological signatures distinct from stan- 
dard weakly interacting dark matter. A modification of 
bubble baryogenesis can achieve simultaneous generation 
of baryons and dark matter by extending the symmetry 
structure to U(1)b x U(1)dm- 
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